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Abstract 

ObsCTvations  regarding  the  use  of  advanced  computational  fluid  dynamics  (CFD)  analysis 
sensitivity  analysis  (S  A),  and  design  codes  in  gradient-based  multidisciplinary  design  optimization 
(MDO)  reflect  our  pox^eption  of  the  interactions  required  of  CFD  and  our  experience  in  recent 
aerodynamic  design  optimization  studies  using  CFD.  Sample  results  from  these  latte-  studies  are 
summarized  for  conventional  optimization  (analysis-SA  codes)  and  simultaneous  analysis  and 
design  optinuzation  (design  code)  using  both  Euler  and  Navier-Stokes  flow  {q)proximations.  The 
amount  of  computational  resources  required  for  aerodynamic  design  using  CFD  via  analysis- 
SA  codes  is  greater  than  that  required  for  design  codes.  Thus,  an  MDO  fonnulation  that 
utilizes  the  more  efficient  design  codes  where  possible  is  desired.  However,  in  the  aerovehicle 
MDO  problem,  the  various  disciplines  that  are  involved  have  different  design  points  in  the  flight 
envelope;  therefore,  CFD  analysis-SA  codes  are  required  at  the  aerodynamic  “off  design”  points. 
The  suggested  MDO  formulaticm  is  a  hylxid  multilevel  optimization  procedure  that  consists 
of  both  multipoint  CFD  analysis-SA  codes  and  multipoint  CFD  design  codes  that  pafoim 
suboptimizations. 
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1  Introduct<H‘y  Remarks 

The  focus  of  this  work  is  on  only  those  techniques  that  are  applicable  to  advanced  (high-fidelity) 
computational  fluid  dynamics  (CFD)  and  that  are  extendable  to  multidisciplinary  design  optimization 
(MDO)  for  realistic  three-dimensional  (3-D)  aerovehicles.  Table  1  briefly  summarizes  this  research  and 
gives  the  chief  advantages  and  disadvantages  of  each  technique. 


Table  1.  Summary  of  Techniques  Used  in  Present  Focus 


TECHNIQUE 

ADVANTAGE 

DISADVANTAGE 

Advanced  CFD 

Required  physical  fidelity 

Nonlinear  and  high  costs 

Multiblock  (or  unstructured) 

Required  geometric  complexity 

Higher  costs 

Multigrid  acceleration 

Efficient  algcHithm 

Tedious  code 

Gradient-based  optimization 

gradient  informaticxi  in  math 
models 

Local  vs.  global  minima 

Automatic  diff^entiation 

Accurate  and  robust 

Efficiency  questions 

Incremental  iterative  method 

Ctxisistent,  efficient,  and 
versatile  SA  algtmthm 

Not  automatic 

Application  of  these  techniques  to  advanced  CFD  codes  has  been  proposed  and  discussed  in  [l]-[5].  In 
[1]  and  [2],  the  incr^ental  iterative  method  (HM)  for  calculation  of  sensitivity  derivatives  (SD’s)  is 
discussed.  In  [3],  the  {application  of  automatic  differentiation  (AD)  to  obtain  SD’s  from  a  3-D  thin-layer 
NaviCT-Stokes  code  is  demonstrated.  References  [4]  and  [5]  are  rec«it  summaries  of  the  studies  of  this 
work;  the  combination  of  AD  and  EM  to  efficiently  obtain  consistent  discrete  SD’s  frtMn  a  two-dimensional 
(2-D)  thin-layer  Naviar-St(*es  code  is  demonstrated  in  [6].  The  AD  tool  ADIFOR  (AD  of  Fortran)  of  [7] 
and  [8]  has  been  used  throughout  this  work.  Symbols  and  acnmyms  are  defined  as  introduced  in  the  text 
The  equations  summarized  in  this  paragr£q)h  can  be  found  in  greater  detail  in  [l]-[6].  The  conservatitm 
laws  of  comivessible  fluid  flow,  R,  and  aerodynamic  functions,  F,  of  interest  can  be  expressed  as 

(I)  i2(C?(6),  X(6),  6)  =  0  {N onlinear  state  equation) 

and 


(2)  F  =  F(Q(b),X(b),b)  (Aerodynamic  output  function) 

where  (2  is  the  vector  of  state  (field)  variables,  X  is  the  vector  of  computational  grid  coordinates,  and 
b  is  the  vector  of  design  variables.  Direct  differentiation  of  Eqs.  (1)  and  (2)  with  respect  to  the  design 
variables  yields 


(3) 

and 


R'  = 


dR 

db 


(4) 


dF 

dQ 


dF 

db 


where  Q'  =  ^  and  X'  =  ^  are  the  flow  and  grid  SD’s  with  respect  to  the  design  variables.  The 
introduction  of  an  adjoint  variable  (vector)  A  associated  with  F  gives 


(5) 


{Costate  equation) 


1 


and 


dR 


(6)  ‘  dh  '  56 

where  T  denotes  transpose.  The  IIM  solution  forms  for  Eqs.  (1),  (3),  and  (5)  are,  respectively, 


(7) 


dQ 


AQ  =  i?"; 


g"+i  =  Q”  +  A(?  (n  =  l,2,3,...) 


(8) 


-|§A<J'  =  ir-; 


Q/(m+l)  ^  g/(m)  ^  AQ'  (m  =  1,  2,  3, ...) 


whoe  n  and  m  are  iteration  indices  and  |§  denotes  an  “^proximate  operatOT  of  convenience.”  The  first- 
order  SD’s  (r)  are  obtained  firom  either  Eq.  (4)  with  O'  and  X'  or  Eq.  (6)  with  A  and  X'.  Second-order 
SD’s  are  also  discussed  in  [6]  but  not  used  in  the  present  work. 

Calmlation  of  SD’s  via  finite  differences  (FD)  for  large  sets  of  iteratively  solved  nonlinear  equations 
can  be  expensive  and  inaccurate.  A  comparison  of  SD  ratios  for  lift,  drag,  and  pitching  moment  with 
respect  to  a  geometric  design  variable  fm  transonic  turbulent  flow  about  an  airfoil  [9]  is  shown  in  Table  2. 
The  tiftnnminarnr  in  these  SD  ratios  is  the  value  of  the  SD  via  AD  at  relative  residual,  R/ Ri,  and  residual 
derivatives,  (^)/(^)  i.  reductions  of  10~^.  If  the  FD  results  (SDfd)  agreed  with  the  AD  results  (SDad), 


Ihble  2.  Sensitivity  Derivative  Ratio  Comparisons:  A  Turbulent  Viscous  Example 


R 

%  change  of  design  variable,  b 

(^) 

SDad 

10 

1 

0.1 

0.01 

0.001 

/  dR\ 

\  / 1 

{SDad)io-6 

lift 

0.937 

0.926 

0.040 

— 

— 

10“^ 

drag 

4.862 

1.819 

4.383 

— 

— 

10“^ 

1.153 

moment 

1.035 

2.097 

-0.099 

— 

— 

lift 

1.043 

0.973 

0.873 

-0.018 

— 

0.999 

1 

o 

drag 

3.792 

2.585 

7.883 

29.154 

— 

lO"'^ 

1.007 

moment 

1.139 

0.944 

1.129 

0.385 

— 

0,999 

lift 

— 

0,^5  ^ 

1.039 

0.973 

-0.003 

0.999 

10“5 

drag 

— 

U67  : 

-0.159 

6.986 

28.851 

tn 

1 

O 

1.000 

moment 

— 

om  \ 

1.188 

1.269 

0.397 

0.999 

lift 

— 

1.003 

0.992 

1.112 

1.331 

1.000 

10“* 

drag 

— 

1.507 

0.991 

-0.904 

2.804 

10“« 

1.000 

moment 

— 

1.006 

1.003 

1.270 

1.680 

1.000 

all  table  entries  would  be  like  those  in  the  last  column  (i.e.,  essentially  unity).  The  table  shows  that  at  a 
given  convogence  level  for  R,  different  output  functions  may  require  different  design  variable  step  sizes 
(and  vice  versa)  for  the  SDfd-  The  similarly  shaded  areas  show  comparable  SD  accuracies  for  the  FD 
and  AD  tqrproaches. 
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Recent  ovoviews  of  optimization  ^>proaches  for  coupled  systems  (disciplines)  are  given,  for  example, 
in  [10]  and  [1 1];  the  present  single-discipline  optimization  observations  will  be  discussed  in  terms  similar  to 
those  introduced  in  [10].  Ctxiventional  optimization  utilizes  analysis  (and  perh^s  also  sensitivity  analysis 
(SA))  codes  in  conjuncti(Hi  with  an  optimization  code,  as  shown  in  Fig.  1(a);  when  an  analysis  is  an  iterative 
one,  it  must  be  run  to  reasonably  well-converged  solutions  many  times.  In  [10],  this  method  is  called  the 
nested  analysis  and  design  (NAND)  approach;  the  optimizer  and  analysis  codes  are  loosely  coupled.  At 
the  other  extreme  is  the  simultaneous  analysis  and  design  (SAND)  approach  in  which  the  design  and  state 
variables  are  updated  together.  For  an  iterative  (stale)  analysis  code,  then,  the  design  variable  updates  are 
made  within  the  iterative  analysis  loops;  the  optimizer  and  analysis  are  implicitly  coupled  to  produce  an 
iterative  design  code  as  shown  in  Fig.  1(b).  In  reality,  an  entire  range  of  design  procedures  exists  between 
the  two  extremes  of  Figs.  1(a)  and  1(b);  these  procedures  differ  only  in  the  frequency  at  which  the  iterative 
analysis  and  optimization  interact 


Conventional 


I  optimizerH 


giv« 

«is 

izr“ 


Objective  and  constraints 
De^gn  variables 

Loosely  coupled 


Equilibrium  conditions 
State  variables 


HAND 


Simultaneous  Analysis  and  Design 
Optimization 


Simultaneous 

update 


SAND 


(a)  Nested  analysis  and  design  (NAND).  (b)  Simultaneous  analysis  and  design  (SAND). 

Fig.  1.  Approaches  for  an  optimization  with  iterative  (analysis)  state  equation. 


In  the  next  few  sections,  this  paper  addresses  recent  aerodynamic  single-point  design  optimization 
studies  using  both  NAND  and  SAND  approaches,  with  a  comparison  of  computational  requirements; 
proposed  aerodynamic  multipoint  aerodynamic  design  optimization  approaches;  and,  finally,  a  suggested 
MDO  ^proach. 


2  CFD  Single-Pmnt  Design  Optimization  Results 

The  results  discussed  in  this  section  are  for  recent  single-discipline  (CFD)  single-point  design  optimization 
studies;  details  have  been  reported  elsewhere  as  noted  in  the  quoted  references.  Both  NAND  and  SAND 
^proach  results  have  been  obtained,  although  generally  not  for  identical  problems.  Even  though  formal 
optimization  procedures  have  been  used,  these  studies  produce  design  improvements  as  opposed  to  optimum 
solutions,  probably  because  of  the  accepted  convergence  levels  in  the  required  iterative  solutions  for  the 
nonlinear  CFD  equations. 

2.1  Conventional  Optimization:  NAND  with  SA 

The  conventional  optimization  study  results  discussed  here  are  for  aerodynamic  sh^  optimization  using 
the  NAND  ^proach  with  SA  based  on  analytical  SD’s.  The  High-Speed  Civil  Transport  (HSCl)  24E 
design  improvement  studies  are  based  on  a  3-D  marching  (supersonic)  Euler  CFD  code,  and  results  are 
reported  in  [12].  The  transonic  turbulent  airfoil  study  results  are  based  on  a  2-D  thin-layer  Navier-Stokes 
CFD  code,  and  results  are  reported  in  [9]. 

2.1.1  HSCT  Design  Improvement.  Initial  results  for  aerodynamic  shape  optimization  studies  that  in¬ 
vestigate  the  feasibility  of  using  a  3-D  supersonic  Euler  code  with  an  efficient  SA  capability  are  given  in 
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[12]  fOT  a  Mach  2.4  HSCT  wing-body  configuration.  A  comparison  of  typical  nongeomertc  SD  s  of  ^ 
f(Mce  and  moment  coefficients  with  respect  to  Mach  number,  angle  of  attack,  and  yaw  angle  is  given  m  [  ]. 
The  nM  and  FD  SD’s  agree  to  four  significant  digits;  the  IIM  results  are  computationally  less  expensive 
to  obtain.  TVpical  aerodynamic  SD’s  of  the  force  and  moment  coefficients  with  respect  to  wing  geoinetry 
parameters  obtained  with  the  IIM  are  also  accurate  and  computationaUy  less  expensive  than  those  obtained 
by  efficient  FD.  A  comparison  of  these  geometric  SD’s  with  respect  to  wing-section  thickness,  twist,  and 
cambCT,  as  well  as  wing  planform  and  flap  deflections  is  given  in  [12].  Again,  the  agreement  between  the 
nM  and  FD  results  is  very  good,  and  the  IIM  results  are  computationally  less  expensive  to  obtain. 

The  cOTiputational  flowchart  for  the  aerodynamic  shj^  optimization  studies  is  shown  m  Fig.  2, 
where  the  (outer)  shape  design  iteration  loop  is  shown.  In  these  studies,  extensive  use  was  made  of 
solution  restart  files;  these  inner  loops  are  not  shown  in  Fig.  2.  The  shape  design  loop  starts  at  the  u^r 
left  with  the  autcmated  surface-shying  and  volume-grid  generatiCHi  codes,  which  are  discussed  m  [1  ] 
and  [15].  Tliese  codes  are  differentiated  with  ADIFOR  ([7]  and  [8])  to  provide  the  grid  SD’s  (X )  with 
respect  to  jrairoximately  100  (wing)  geometric  design  variables,  as  discussed  in  the  ^pendix  of  [12]. 
Both  the  grid  (X)  and  its  SD’s  (X')  are  required  because  the  geometric  (sh^)  design  variables  determine 
the  vehicle  surface  and  its  body-fitted  computational  grid.  The  marching  Euler  code  is  differentiated  by 
hand  (|§,  If,  and  and  likewise  for  ouqiut  functions  F)  to  construct  the  flow  derivative  code.  The 
AutomS  Design  Synthesis  (ADS)  program  [16]  is  used  for  the  present  constrained  optimizarion  results; 
the  sequential  quadratic  programming  strategy,  the  modified  method  of  feasible  directions  optimize,  m 
the  Golden  Section  line  search  options  have  been  selected.  Evaluation  of  both  function  and  firsttorder 
derivatives  (SD’s)  is  given  to  the  ADS  code.  Because  the  SD  via  the  IIM  are  essentially  analytical 
derivatives,  this  combination  of  methods  in  ADS  gives  the  most  consistent  optimization  results.  However, 
many  function  evaluations  are  required  by  the  selected  search  procedure.  Increments  in  the  sh^  design 
variables,  draioted  “Del  geom  input”  in  Fig.  2,  are  returned  to  the  surface  shaping  code  to  start  the  next 

design  iteration. 


Fig.  2.  Flowchart  for  aerodynamic  shape  optimization  feasibility  studies. 

Sample  results  from  two  of  the  eight  optimization  studies  rqported  in  [12]  are  discussed  here.  The 
HSCT  24E  filleted  wing-body  configuration  generated  by  NASA  Langley  Research  Center  is  the  b^elme 
shye;  the  flow  conditions  are  Mach  number  Mo©  =  2.4,  angle  of  attack  =  1®,  and  yaw  angle  /?  -  . 

Conv«gence  of  both  the  nonlinear  iterative  flow  analysis  and  the  linear  iterative  SA  is  to  a  relative 
residual  reduction  of  6  orders  of  magnitude  for  all  required  solutions. 

For  the  wing-section  thickness  design  improvement  study,  initial  and  final  thickness  distributions  are 
shown  in  Fig.  3.  The  15  design  variables  consist  of  5  parameters  each  at  the  wing  root,  bre^,  and 
tip  locations.  The  wing  is  linearly  lofted  frcMn  root  to  break  and  fr(Mn  break  to  tip  to  supply  thickness 
information  at  all  othw  wing  stations.  The  objective  function  is  drag  minimiration,  with  the^wing  r^ 
bending  moment  and  lift  constrained  to  their  baseline  values;  that  is,  minimize  ^  subject  to  <  1.0 
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^  l-O-  The  baseBne  drag  is  decreased  by  about  10  pwcent,  and  both  constraints  are  active.  This 
improvement  is  obtained  in  8  optimization  steps,  which  requires  1 17  fimction  evaluations  and  8  gradient 
evaluations;  the  Cray-2  run  time  is  about  1.2  hours.  For  6  of  the  15  design  variables,  the  side  constraints 
are  active  (vnthin  5  pocent  of  the  specified  bounds,  which  were  arbitrarily  taken  for  the  thickness  variables 
to  be  ±50  percent  of  the  baseline  values).  For  sup^sonic  flow  considerations  alone,  the  wing  would  be 
expected  to  become  thinner,  as  shown  in  Fig.  3. 


Fig.  3.  Wing  section  thickness  distributions  for  HSCT  24E  design  improvement. 

For  the  wing  planform  design  improvement  study,  initial  and  final  planforms  are  shown  in  Fig.  4. 
The  five  design  variables  are  the  wing  root,  break,  and  tip  chord  lengths  and  the  streamwise  locations  of 
both  break  and  section  leading  edges.  The  objective  is  to  minimize  —■§f  -  subject  to  <1.0  and 

^  <  1.0.  The  baseline  lift  is  increased  by  5.5  percent,  and  the  drag  constraint  is  violated  by  3.8  percent 
Neither  the  wing  bending-moment  constraint  nor  any  of  the  design-variable  side  constraints  are  active  or 
violated.  For  supersonic  flow  considerations  alone,  the  wing  tip  should  be  swept  more  than  in  the  baseline 
HSCT  24E;  Fig.  4  shows  that  the  optimization  procedure  produces  that  result  At  a  Mach  number  of  2.4, 


Fig.  4.  Wing  planform  shape  for  HSCT  24E  design  improvement. 

the  Mach  angle  {pi)  is  24.6°.  The  angle  subtended  by  the  wing-tip  leading  edge  from  the  root  leading  edge 
is  25.9  for  the  baseline  HSCT  24E  and  23.8°  for  the  final  optimized  planform.  That  is,  the  planform 
optimized  for  only  supersonic  flow  lies  behind  the  Mach  cone. 

2.1,2  Ttirbulent  Transonic  Airfoil  Improvement.  Initial  results  for  aerodynamic  shape  optimization 
studies  that  investigate  the  feasibility  of  using  a  2-D  Navier-Stokes  code  with  an  efficient  SA  capability 
are  given  in  [9]  for  a  turbulent  transonic  airfoil.  Comparisons  of  typical  SD’s  of  force  and  moment 
coefficients  with  respect  to  both  flow  and  geometric  variables  are  given  in  [6]  for  several  methods  of 
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Cflirniatinn  The  agreement  between  the  EM  via  AD  and  FD  is  good;  the  computational  costs  for  the  IIM 
using  AD  are  comparable  to  or  less  than  those  for  the  FD. 

The  flowchart  for  this  airfoil  design  improvement  process  is  essentially  that  shown  in  Fig.  2;  a  more 
detailed  flowchart  with  the  restart  loops  indicated  is  given  in  [9].  The  automated  geometry  and  grid  code 
TBGG  (two-boundary  grid  generation)  [17]  is  used  for  algebraically  generated  C  grids  around  airfoils. 
The  airfoil  upp^  and  lower  surfaces  are  each  represented  by  a  linear  combination  of  four  orthonormalized 
polynomials  [18];  their  amplitude  (or  weights)  are  the  sht^  design  variables.  The  entire  TBGG  code  is 
differentiated  with  ADIFOR  ([7]  and  [8])  to  obtain  the  grid  SD’s  (X')  with  respect  to  these  eight  design 
variables.  The  CFD  code  ANSERS  (algcwithm  for  the  Navier-Stokes  equations  based  on  a  Riemann 
solver)  [19]  is  a  2-D  Navier-Stokes  code  based  on  an  upwind  cell-centered  finite-volume  formulation 
and  simulates  turbulence  wifli  the  Baldwin-Lomax  [20]  algebraic  turbulence  model.  This  code  is  also 
differentiated  with  ADIFOR;  however,  the  differentiation  is  done  in  parts  so  that  an  efficient  IIM  can  be 
constructed  for  the  SD’s  [6].  As  in  the  previous  example,  the  ADS  program  [16]  is  used  for  the  present 
ron«;ti^in<»/i  optimization  results.  The  sequential  quadratic  programming  strategy,  the  modified  method  of 
feasible  directions  optimize,  and  the  Golden  Section  line  search  optirxis  are  selected;  evaluation  of  both 
the  function  and  first-order  derivatives  is  provided  to  the  ADS  code. 

The  design  optimization  problem  summarized  here  (and  fully  discussed  in  [9])  is  the  maximization 
of  the  lift-to-drag  ratio  of  an  airfoil  in  a  turbulent  transonic  flow.  The  incident  flow  is  at  a  Mach  numba- 
of  0.8,  a  chord  Reynolds  numba  of  5  million,  and  an  angle  of  attack  of  1®.  Solutions  are  obtained  on  a 
C-^  mesh  of  257  (circumferential  direction)  x  67  (nwmal  direction)  points.  In  order  to  maintain  the 
convexity  of  the  airfoU  around  the  leading  edge,  10  constraints  of  positive  curvature  are  imposed  at  equally 
spaced  grid  points  on  both  surfaces  for  the  first  5  percent  chord. 

The  initial  solution  (Q)  for  the  NACA  0012  airfoil  is  converged  to  a  relative  analysis  residual  reduction 
of  10  orders  of  magnitude,  whereas  for  the  SD  (Q')  the  relative  derivative  residuals  are  converged  6  orders. 
These  solutitms  ware  used  as  initial  values  for  design  optimization  iterations.  An  optimal  solution  was 
obtained  after  three  resubmittals  to  ADS;  these  resubmissitms  ^sure  that  the  optimization  iterations  have 
not  terminated  iwerrraturely  and  are  denoted  by  the  black  symbols  on  the  objective  function  history  plot 
of  Fig.  5.  The  objective  function  Ct/Cd  has  been  increased  to  ^proximately  10  times  its  initial  value  in 
58  design  iterations.  This  optimization  requires  the  same  amount  of  central  processing  unit  (CPU)  time 
as  required  for  tqrproximately  36  well-cwivaged  single  aerodynamic  analyses.  The  strrmg  shocks  present 
in  the  initial  design  have  been  eliminated  from  the  final  design,  as  can  be  seen  from  the  surface  pressure 
distribution  plots  (Cp)  in  Fig.  5  and  the  flow-field  Mach-numbw  contour  plots  shown  in  Fig.  6.  The  airfoil 
has  been  thinned  and  cambwed  considerably,  as  would  be  expected  in  order  to  maximize  the  lift-to-drag 
ratio  at  transcxiic  flow  conditions. 


Fig.  5.  Objective  function  and  surface  pressure  distributions  for  turbulent  transonic  airfoil  improvement. 

In  the  results  above,  the  relative  residual  reductions  for  both  Q  and  Q'  were  set  to  4  orders  in  the 
design  optimization;  addirinnal  solution  iterations  of  Eq.  (1)  were  required  to  improve  the  solution  quality 
of  the  final  design.  However,  a  limited  study  of  several  cases  to  assess  the  effects  of  both  flow  analysis 
and  SA  convwgence  tolerances  on  the  performance  of  the  aerodynamic  design  optimization  is  given  in  [9]. 
Different  cf>mhinatir>ns  of  flow  and  SA  tolerances  (differing  by  factors  10*^*)  result  in  different  achieved 
objective  functions  (differing  by  50  percent  to  75  percent);  these  diffCTent  combinations  require  different 
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total  CPU  times  and  produce  different  airfoils.  Such  different  results  are  not  unexpected  in  this  gradient- 
based  optimization  procedure,  because  the  (partially  converged)  values  of  aerodynamic  functions  and 
d«ivatives  govern  the  (^proximate)  objective  function  surface  and  the  (tpproximate)  search  direction  mi 
that  surface.  The  different  airfoils  that  result  may  also  correspond  to  different  local  optima  on  the  true 
(well-CMivMged)  objective  function  surface. 


(a)  Initial.  (b)  Final. 

Fig.  6.  Flow-field  Mach  number  contour  plots  for  turbulent  transonic  airfoil  improvement, 
in  Simultaneous  Analysis  and  Optimization:  SAND 

Several  simultaneous  analysis/optimization  demonstration  study  results  are  discussed  here  for  aerodynamic 
sh^  optimization  using  a  SAND  ^proach  with  an  SA,  based  on  analytical  SD’s  or  adjoints.  The 
present  procedure,  using  CFD,  has  been  called  SAADO  (simultaneous  aerodynamic  analysis  and  design 
optimization).  This  procedure  incMporates  design  optimization  within  CFD  analysis  to  achieve  a  convMged 
flow  solution  and  an  optimal  design  at  the  same  time.  For  advanced  3-D  CFD  codes,  an  iterative  solution 
of  the  linearized  ^)proximations  to  the  nonlinear  flow  equations  is  required  because  of  the  large  matrirpc 
involved.  Design  optimization  is  also  iterative;  the  SAADO  procedure  intoacts  these  two  iterations 
(before  either  is  conv«ged),  allowing  a  simultaneous  relaxation  of  both  the  flow-field  solution  and  the 
design  optimization.  Ovoall  computational  efficiency  is  achieved  because  expensive  iterative  solutions 
for  nonoptimal  design  parameters  are  not  converged  (i.e.,  obtained).  Initial  demonstration  studies  are 
for  quasi-one-dimensiMial  (1-D)  nozzle  flows  described  by  an  Euler  equation;  the  results  are  reported  in 
[21]  and  [22].  The  tiansMiic  turbulent  airfoil  results  are  based  on  a  2-D  thin-layer  Navier-Stokes  CFD 
tqilHOximation  and  the  results  are  reported  in  [23]. 


2.2.1  Quasi-1— D  Nozzle  Demonstration.  The  basic  formulation  and  several  variations  of  SAADO 
have  been  derived  and  successfully  implemented  for  design  optimization  of  quasi-l-D  nmyipt  [21,  22]. 
The  feasibility  of  SAADO  is  demonstrated  for  both  supersonic  and  transonic  flow  described  by  a  1-D 
Euler  equation.  OptimizatiMi  results  for  a  supersonic  nozzle  design  are  shown  in  Fig.  7(a).  Standard 
NAND  optimization  results  with  ADS  software  [16],  for  both  one-sided  FD  and  quasi-analytical  (QA)  [2] 
differentiation  for  SD’s  are  compared  with  the  SAADO  results.  The  SAADO  procedure  produces  a  better 
objective  minimization  (which,  in  this  case,  should  tend  to  zero)  in  20  design  cycles  with  fewer  equivalent 
Newton-Raphson  (NR)  iterations  and  in  less  computer  time  than  either  of  the  standard  optimizations.  For 
the  transMiic  nozzle,  the  residual  (R)  history  plots  in  Fig.  7(b).  show  that  SAADO  requires  approximately 
the  same  number  of  iterations  to  reach  the  optimal  solution  as  is  required  for  a  single  flow  analysis.  The 
osciUating  nature  of  the  SAADO  residual  history  indicates  that  the  incorporation  of  design  changes  alters 
the  normal  convergence  pattern  for  flow  analysis;  these  oscillations  cease  when  the  design  changes  cease. 
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In  both  nozzle  designs,  the  objective  is  to  match  a  prescribed  velocity  distribution  throughout  the  length 
of  the  nozzle.  Ihe  results  are  shown  in  Fig.  8. 


Optfmizer 

ADS 

ADS 

mm 

SDcalculalion 

FD 

QA 

Objective 

8.8E-4 

8.9E-4 

Equiv.  NR  iter. 

174 

102 

Cray-2  (sec) 

2.19 

1.29 

{a)  Supersonic. 


Fig.  7.  Convergence  results  for  SAADO  quasi-1^  nozzle  flow  feasibility  demonstration. 


%  of  duct  length 


Fig.  8.  Nozzle  velocity  distributions  for  supersonic  and  transonic  flow  SAADO  feasibility  demonstration. 

The  primary  objective  of  the  SAADO  formulation  is  a  significant  reduction  in  the  CPU  time;  the 
CPU  time  required  for  the  repeated  analyses  in  standard  NAND  optimization  is  impractical.  The  potential 
for  SAADO  is  clearly  demonstrated  in  these  sample  1-D  nozzle  design  optimization  problems.  For  the 
supersonic  flow  case,  the  better  solution  is  obtained  by  SAADO  in  53  pwcent  of  the  CPU  time  required  by 
the  NAND  (QA  for  SD)  ^proach;  this  solution  time  corresponds  to  the  time  required  for  approximately 
9  equivalent  analyses.  For  the  transonic  flow  case,  the  SAADO  approach  obtains  the  design  in  the  time 
required  for  1  or  2  analyses. 


2.2,2  Tiirbulent  IVaiisonic  Airfoil  Demonstration.  The  initial  results  for  aerodynamic  shape  optimiza¬ 
tion  with  a  2-D  SAADO  procedure  include  both  Euler  and  thin-layer  Navier-Stokes  CFD  tq)proximations. 
Results  for  the  latter  are  reported  in  [23]  and  are  briefly  discussed  here.  The  physical  problem  is  de¬ 
scribed  in  section  2.1.2;  howevCT,  airfoil  thickness  constraints  have  been  added  to  account  for  a  structural 
wing  box.  This  2-D  SAADO  procedure  is  constructed  ffCHn  the  same  grid-generation  (TBGG  [17])  and 
flow  (ANSERS  [19])  codes  briefly  discussed  in  section  2.1.2.  Both  hand  differentiation  (Euler,  adjoint) 
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and  ADIFOR  (Navier-Stokes,  direct  differentiation)  have  been  used  to  generate  the  required  sensitivity 
equations  for  the  SAADO  versions  that  have  been  tested. 

Presently,  the  2-D  form  of  the  SAADO  procedure  for  a  thin-layer  Navier-Stokes  code  consi^rs 
both  airfoil-shape  and  flow-field  variables  as  independent  (design)  variables  and  treats  the  flow  equations 
as  equality  constraints.  This  large  number  of  variables  and  constraints  is  reduced  by  solving  a  set  of 
inexact  sensitivity  equations;  these  equations  and  the  flow  analysis  are  solved  together  in  IIM  form.  The 
design  optimization  problem  is  to  maximize  the  lift-to-drag  ratio  CtjCd  of  an  airfoil  in  transonic  turbulent 
flow,  subject  to  25  geometric  constraints:  20  on  the  surface  curvature  imposed  from  the  leading  edge  to 
5  percent  chord  and  5  on  the  airfoil  thickness  imposed  from  20  to  60  percent  chord.  Airfoil  upper  and 
lower  profiles  are  each  represented  by  a  linear  combination  of  four  orthonormalized  polynomials,  as  in 
[18].  The  initial  values  of  these  eight  polynomial  weighting  coefficients  (the  sh^  design  variables)  are 
those  for  the  NACA  0012  airfoil. 

Sample  results  for  the  SAADO  approach  to  improving  the  lift-to-drag  ratio  at  a  Mach  number  of 
0.8,  an  angle  of  attack  of  1°,  and  a  chord  Reynolds  number  of  5  million,  starting  from  the  NACA  0012 
airfoil,  are  obtained  on  a  C-type  mesh  of  127  x  33  points.  Figure  9  shows  the  objective  function  (versus 


Final 

Initial 


Fig.  9.  Objective  function  and  surface  pressure  distributions 
for  turbulent  transonic  airfoil  SAADO  demonstration. 

number  of  SAADO  cycles)  history  and  surface  pressure  (coefficient,  -Q?)  distributions  on  the  initial  and 
final  airfoil  surfaces.  In  18  SAADO  cycles  (design  variable  updates),  the  flow  analysis  convergence  level 
was  reduced  firom  10~^  to  10“^,  which  produces  an  airfoil  with  ^proximately  500  percent  improvement  in 
the  objective  function.  Two  thickness  constraints,  (at  20  and  60  percent  chord)  were  active.  These  results 
show  that  the  shock  wave  on  the  lower  surface  of  the  airfoil  has  been  eliminated;  the  shock  wave  on  the 
upper  surface  has  been  weakened  considerably.  These  different  shock  wave  patterns  can  also  be  seen  in 
the  initial  and  final  flow-field  Mach  number  contour  plots,  which  are  shown  in  Fig.  10. 

The  computational  time  required  for  this  SAADO  case  is  approximately  the  same  as  22  flow  analyses 
converged  to  a  relative  residual  error  of  10“^  on  this  127  x  33  grid.  For  the  inviscid  (Euler)  approximaticHi, 
the  relative  SAADO  time  has  been  reduced  further.  The  set  of  eight  inexact  aerodynamic  sensitivity 
equations  is  replaced  by  three  previously  obtained,  hand-differentiated  adjoint  equations  reducing  the 
computational  time  by  ^proximaiely  one-third.  For  a  single  adjoint  equation,  which  corresponds  to 
the  single  ouqiut  function  C^/C^,  the  estimated  reduction  in  computational  time  would  be  by  a  factor  of 
^proximately  3.  This  reduction  is  not  currently  feasible  for  the  Navier-Stokes  code  because  the  ADIFOR 
tool  does  not  provide  a  computationally  efficient  adjoint  equation  code.  Nevertheless,  the  feasibility  for 
implementation  and  use  of  the  SAADO  procedure  on  a  2--D  thin-layer  Navier-Stokes  code  in  transcmic 
turbulent  airfoil  design  improvement  has  been  demonstrated;  ADIFOR  was  used  to  generate  the  derivative 
code. 

23  Comments  on  NAND  versus  SAND  Optimization. 

A  number  of  previous  formulaticms  that  differ  frcnn  SAADO  have  also  involved  efficient  SAND  method¬ 
ologies  for  aerodynamic  design  optimization.  The  works  of  Rizk  [24],  Campbell  [25],  Drela  [26],  and 
Young  et  al.  [27],  for  example,  either  incorporated  simple  design  modification  rules  in  the  aerodynamic 
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{a)  Initial.  (b)  Final. 

Fig.  10.  Flow-field  Mach  number  contour  plots  for  turbulent  transonic  airfoil  SAADO  demonstration. 
analysis  iterations  or  extended  the  Newton-like  methods  to  include  the  iteration  equation  that  guides  the 
change  of  design  variables.  Ghattas  and  Orozco  [28]  and  Ta’asan  et  al.  [29]  also  developed  procedures 
for  SAND.  Ghattas  and  Orozco  derived  an  equation  that  relates  the  design  changes  to  the  changes  in  the 
flow  solution,  based  on  the  sparsity  of  the  Hessian  mauix  in  the  design  optimization  problem  formulatirai. 
Ta’asan  et  al.  integrated  the  aerodynamic  design  optimization  within  the  multigrid  method  for  aerodynamic 
analysis.  The  error  reduction  in  each  grid  improves  not  only  the  quality  of  the  flow  solution  but  also  the 
design  improvement  at  that  scale.  Results  from  other  SAND  procedures  are  given  in  several  p^)ers  also 
contained  in  this  volume. 

The  point  here  is  that  SAND  procedures  produce  aerodynamic  designs  much  more  efficiently  than  the 
NAND  procedures.  Solution  times  for  the  SAND  procedures  are  measured  in  one  to  a  few  single-analysis 
times,  wh^eas  those  for  even  v^  efficient  NAND  procedures  are  measured  in  many  equivalent  single¬ 
analysis  times.  This  diffi^nce  was  shown  for  the  present  nozzle  example  discussed  in  section  2.2.1,  where 
the  SAND  solutions  varied  from  tq)proximateiy  2  to  9  equivalent  single-analysis  times  and  the  NAND 
solutions  woe  another  factor  of  2  to  3  times  larger.  For  the  turbulent  transonic  airfoil  examples  discussed 
in  sections  2.12  (NAND)  and  2.2.2  (SAND),  this  difference  is  not  as  clear  because  structural  wing-box 
thickness  constraints  w^  imposed  only  in  the  latter  example  and  the  computational  grids  were  different.  In 
addition,  the  SAADO  was  implemented  with  ADIFOR  (fcHward-mode)  direct  differentiation  rather  than  the 
(revQ^-mode)  more  efficient  adjoint  formulatirai.  Nevertheless,  the  efficient  NAND  procedure  required 
^proximately  36  equivalent  single-analysis  solution  times;  the  SAND  procedure  required  22  single-analysis 
times.  Based  on  Euler  solution  comparisons,  the  estimated  time  for  a  single-adjoint  SAND  procedure  for 
the  Navier-Stokes  equations  would  be  equivalent  to  7  single-analysis  times. 


3  CFD  Multipoint  Design  Optimization  Plans 

Multipoint  design  for  any  single-discipline  optimization  generally  exhibits  several  facets  of  an  MDO 
formulatim.  Features  of  the  solutions  at  the  multiple  design  points  may  appear  as  constraints  or  objectives 
(or  parts  of  them)  for  the  multipoint  problem  and  must  be  jointly  optimized  and/or  coordinated  (which  also 
may  be  posed  as  an  optimization)  in  the  space  of  multipoint  design  variables.  This  set  of  multipoint  design 
variables  may  also  contain  subsets  of  variables  for  the  different  single-point  problems.  Thus,  solutions  at 
the  multiple  points  can  be  viewed  and  treated  as  different  disciplines.  When  some  or  each  of  the  multiple 
design  points  involve  a  local  design  or  determine  a  subset  of  the  multipoint  design  variables,  then  multilevel 
optimization  arises  naturally.  Arbitrary  decompositions  of  the  multipoint  design  problem  may  be  based 
on  lower-level  partitionings  guided  by  oth^  considerations  such  as  problem  size,  computer  resources  and 
codes,  domain  decomposition,  etc.;  such  decompositions  lead  to  partition  interface  conditions  which  are 
mismatched  and  must  be  rectified  by  an  upper-level  coordination. 
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Herein,  coordination  tasks  are  viewed  as  those  involved  in  establishing  multipoint  (or  multidisciplinary) 
feasibility;  that  is,  upper-level  omstraints  which  establish  compatibility  of  shared  design  variables,  target 
values,  and  partition  int^ace  conditions  for  the  multiple  single-point  analyses  ca-  sublevel  optimizations. 
The  numba*  of  levels  of  optimizatkm  or  cocxrdination  provide  another  method  for  viewing  multipoint 
design  processes.  A  detailed  description  of  design  vari^les  and  objective  and  constraint  functions  at 
each  optimization  level  depends,  of  course,  on  the  particular  decomposition  of  a  given  multipoint  design 
problem.  Our  intention  here  is  to  show  the  multipoint  nature  of  the  MDO  problem  and  suggest  the  need 
for  efficient  discipline  codes  for  both  analysis  and  design  in  order  to  solve  the  MDO  problem  efficiently, 
without  specification  of  the  decomposition  details.  Formal  SA  may  or  may  not  be  implem^ted  along 
with  the  iterative  analyses;  it  is  not  indicated  in  the  diagrams  shown  in  this  section.  In  the  following 
subsections,  gen^  ideas  related  to  planned  CFD  multipoint  design  optimization  studies  are  used  to 
illustrate  the  points  to  be  made. 


3.1  Single-Level  Optimization 

In  the  single-level  processes,  both  optimization  and  coordination  (OC)  tasks  are  performed  at  the  same 
level.  The  conventional,  or  NAND  procedure,  places  these  OC  tasks  above  or  in  an  outer  iterative  loop 
around  the  multipoint  analyses,  as  shown  in  Fig.  11(a).  The  OC  tasks  provide  design  variables  as  iiq)ut  to 
the  multipoint  analyses,  which  upon  inner  iterative  convergence  provide  the  output  functions  (and  p»h^ 
the  SD  informaticm)  required  by  the  OC  tasks.  A  design  variable  update  then  initiates  the  next  outer  locq) 
iteration.  This  procedure  requires  well-cmiverged  analyses  at  each  of  the  multiple  points  for  every  change 
of  design  variables  or  iteration  step  in  the  outer  (design)  loop. 

The  single-level  simultaneous,  or  SAND  procedure,  can  be  viewed  either  as  incorporating  the  OC 
tasks  within  the  multiple  point  iterative  analysis  loops  or  as  running  a  modified  NAND  procedure  with 
nonconverged  itmo'  or  analysis  loops.  Both  design  and  state  variables  are  simultaneously  updated,  and  the 
multiple  point  analyses  are  not  converged  until  the  desired  design  is  obtained.  As  indicated  in  Fig.  11(b), 
both  state  and  adjoint  (costate  or  derivative)  equations  are  iterated  togeth^,  which  simultaneously  improves 
all  multiple-point  flow  residuals  and  OC  task  objectives. 


MULTIPOINT  ITERATIVE 
DESIGN 


State  end  State  end  State  and 
adjoint  adjoint  adjoint 
1 - 1 - 1 


SAA/D 


(aJNA/^D. 


(b)  SAND. 


Fig.  11.  Multipoint,  single-level  optimization  with  iteratively  solved  state  equations. 


3.2  Multilevel  Optimization 

In  the  multilevel  processes,  optimization  and  coordination  tasks  need  not  be  performed  at  the  same  level; 
suboptimizations  may  be  done  at  a  lower  or  inner  iteration  level.  Subsets  of  design  variables  that  have 
been  essentially  det^mined  at  one  design  point  may  not  be  needed  or  relevant  at  other  design  points,  or 
they  may  only  be  required  to  be  not  incwnpatible.  For  example,  many  design  variables  for  the  high-lift 
system  are  detmnined  at  takeoff  and  landing  conditions  but  are  not  relevant  at  cruise  conditions,  where 
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efficient  flight  conditions  drive  the  detailed  aerodynamic  shape  design.  Design  variables  which  appear  in 
mcHe  than  one  subq)timization  must  be  subject  to  coordination,  generally  at  a  higher  level.  When  more 
than  one  level  of  optimization  exists,  the  procedure  at  each  level  is  specified.  A  NAND-NAND  procedure 
that  involves  CFD  appears  too  expensive  to  be  considered. 

A  multipoint,  multilevel  SAND-NAND  procedure  is  pictured  in  Fig.  12  for  three  sublevel  design 
points.  The  coordination  task  is  assumed  to  be  in  the  top-level  multipoint  SAND  optimizer.  The  iterative 
analysis  at  each  multiple  design  point  is  the  inn^most  of  three  nested  loops  and  may  be  called  many 
times  p^  optimization  cycle,  particularly  if  SA  is  not  used  to  provide  SD  informatics.  Repeated  use  of 
the  iterative  analysis  has  a  dire  effect  when  that  iterative  analysis  is  expensive  CFD.  Design  variabtes 
at  the  top-level  multipoint  SAND  optimizer  genaally  include  target  values  or  paramet^  in  the  sublevel 
optimization  objective  functions.  A^odynamic  shape  optimization  for  the  configuration  or  parts  of  it  at  low 
subsonic,  transraiic,  and  supersonic  flow  conditions  might  be  the  three  sublevel  design  problms  and  the 
resulting  sh^s  must  be  made  compatible  by  the  upper-level  coordination  at  convergence  in  the  outer  loop. 


Fig.  12.  Multipoint,  multilevel  SAND-NAND  optimization  with  iteratively  solved  state  equations. 


A  multipoint,  multilevel  SAND-SAND  prcxedure  is  pictured  in  Fig.  13  for  the  same  three  sublevel 
design  points  as  depicted  in  Fig.  12.  Again,  the  coordination  task  is  assumed  to  be  in  the  tc^-level 
multipoint  SAND  optimization,  where  design  variables  include  target  values  or  parameters  for  the  sublevel 


Fig.  13.  Multipoint,  multilevel  SAND-SAND  optimization  with  iteratively  solved  state  equations. 
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objective  functions.  Use  of  SAND  pnx:edures  for  the  sublevel  iterative  designs  should  be  computationally 
more  efficient  than  the  corresponding  NAND  procedures  with  SA  at  the  sublevel 

33  Hybrid  Optimization 

A  combination  of  the  single-level  and  multilevel  processes  will  be  called  the  hybrid  optimization  process. 
In  Fig.  14,  the  multilevel  part  is  pictured  as  a  SAND-SAND  procedure,  whaeas  the  single-level  part  is 
shown  as  a  NAND  procedure  which  may  or  may  not  include  SA.  Again,  the  coordination  is  assumed  to  be 
at  the  tq)-level  multipoint  optimization,  or  in  the  outer  loop  shown  in  Fig.  14.  In  a  multipoint  aerodynamic 
wing  design  that  utilizes  CFD  techniques,  for  example,  the  three  sublevel  iterative  design  points  (at  the  left) 
may  be  viewed  as  conesponding  to  takeoff  and  landing  (low  speed),  transMiic  cruise,  and  supCTSonic  cruise, 
as  previously  discussed.  The  iterative  analyses  (at  the  right)  may  correspond  to  other  flow  conditions  (i.e., 
aerodynamic  “off-design”  points)  where,  for  example,  wing  root  bending  moment,  shock  strength,  lift,  or 
drag  (or  other  functions)  need  to  be  assessed  at  each  multipoint  design  step  for  constraint  evaluation  or 
other  considerations.  Often,  such  constraints  arise  frOTti  otha-  discipline  considerations  and  are,  therefore, 
required  in  a  realistic  single-discipline  optimization.  Updated  input  at  these  off-design  analysis  points 
would  be  detamined  by  the  coordination  at  the  upper  level  multipoint  design  step. 


Fig.  14.  Multipoint,  multilevel  hybrid  optimization  with  iteratively  solved  state  equations. 


4  Suggested  MDO  Procedure 

A  given  discipline  multipoint  design,  however,  differs  from  that  shown  in  Fig.  14  because  many  of  the 
required  simultaneous  analyses  wUl  be  from  other  disciplines  and  will  generaUy  occur  at  the  off-design 
points  of  those  disciplines.  For  example,  in  the  multipoint  aerodynamic  wing  design  pictured  in  Fig.  15, 
structural  analyses  will  be  required  at  each  sublevel  aerodynamic  design  point  for  a  flexible  wing  and  also 


Fig.  15.  Multipoint,  multilevel  hybrid  aerodynamic  optimization  with  MDA. 

(possibly)  at  the  aerodynamic  off-design  points.  None  of  these  points  may  be  those  that  correspond  to 
the  structural  design  points.  In  addition,  other  discipline  analyses  may  also  be  required  for  the  multipoint 
aerodynamic  optimizer.  The  same  is  true  for  the  corresponding  structural  multipoint  design;  aerodynamic 
analyses  are  required  at  the  “load-case”  flow  conditions,  as  shown  in  Fig.  16.  Generally,  these  points 
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are  not  the  aerodynamic  (stu^pe)  design  points.  Thus,  the  set  of  multipoint  MDO  design  points  spears 
to  be  a  collection  of  hyteid  multilevel  “discipline-design”  multidisciplinary  analysis  (MDA)  procedures, 
each  with  a  diffe-ent  subset  of  design  points  or  flow  conditions.  A  practicable  implementation  of  these 
procedures  requires  efficient  computatitmal  code  for  both  the  discipline  analysis  and  the  iterative  design. 
For  the  expensive  iterative  analyses  of  advanced  CFD,  the  iterative  design  should  be  a  SAND  procedure 
rather  than  a  NAND  procedure  with  SA. 


Fig.  16.  Multipoint,  multilevel  hybrid  structural  optimization  with  MDA. 

The  suggested  MDO  formulaticm  is  a  coordinated  combination  of  hybrid  multilevel  processes  such 
as  those  pictured  in  Figs.  15  and  16.  The  multipoint  hyNid  multilevel  optimizations  will  be  sublevel 
optimizations  in  an  MDO  coordination  level  as  shown  in  Fig.  17.  This  suggested  MDO  formulation 


Fig.  17.  Suggested  MDO  procedure:  coordinated  hybrid  multilevel  suboptimizations. 

is  much  closer  to  the  currently  practiced  design  procedures  than  an  “all  at  once”  MDO  procedure  as 
described,  for  example,  in  [11].  In  effect,  each  discipline  actively  participates  in  its  sublevel  hybrid  design, 
with  supporting  MDA  and  MDO  coordination  from  the  top  level  or  the  outer  loop.  This  suggested  MDO 
formulation  can  also  be  viewed  as  a  multilevel  coordinated  or  a  collaborative  suboptimization  procedure. 


5  Concluding  Remarks 

Observations  have  been  made  regarding  the  use  of  advanced  computational  fluid  dynamics  (CFD) 
analysis,  sensitivity  analysis  (SA),  and  design  codes  in  gradient-based  multidisciplinary  design  optimization 
(MDO).  These  obsmations  reflect  our  perception  of  the  intoactions  required  of  CFD  and  our  experiaice 
in  recent  aerod5mamic  design  optimization  studies  using  CFD.  Sample  results  from  these  latter  studies 
have  been  giunmari^ftd  for  conventional  optimization  (analysis-SA  codes)  and  simultaneous  analysis  and 
design  optimization  (design  code)  using  both  Euler  and  Naviw-Stokes  flow  ^proximations.  The  amount 
of  computational  resources  required  for  aerodynamic  design  using  CFD  via  analysis-SA  codes  is  shown  to 
be  greater  than  that  requited  via  design  codes.  Thus,  an  MDO  formulation  that  utilizes  the  more  efficient 
design  codes  where  possible  is  desired.  However,  in  the  aerovehicle  MDO  problem,  the  various  disciplines 
that  are  involved  have  different  design  points  in  the  flight  envelope;  therefore,  CFD  analysis-SA  codes  are 
required  at  the  aerodynamic  “off-design”  points.  It  is  concluded  that  the  ^propriate  MDO  formulatitMi 
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include  hybrid  multilevel  optimization  procedures  that  consist  of  both  multipoint  CFD  analysis-SA  codes 

and  multipoint  CFD  design  codes  that  perform  suboptimizations. 
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